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This paper consists of two sections. In the first section, we discuss our fits to the latest HiRes monocular 
spectra. We find that the best fit for the extragalactic component has a spectral index of 7 = —2.38 ± 0.04 with 
a distribution of sources varying with a evolution parameter m — 2.8 ± 0.3. In the second section, we discuss 
preliminary results from a new composition measurement using HiRes monocular data. We find a predominantly 
light spectrum above lO^'^'^ eV. 



1. Fitting the HiRes Spectra 

Fitting the ultra high energy cosmic ray 
(UHECR) spectrum is a common way to attempt 
to gain understanding of the sources of these par- 
ticles. These fits often neglect the galactic com- 
ponent of cosmic rays, which is expected to be 
small at the very highest energies. However, the 
energy at which at which this becomes true is not 
known. 

Rather than assume a negligible galactic com- 
ponent, we employ a toy model in which composi- 
tion determines source. In this model, we use the 
recent composition measurements of the HiRes 
Prototype/Ml A experiment ^ and of HiRes in 
stereo mode These are shown in Figure^ The 
light component in these measurements, assumed 
to be protons, is identified with the extragalactic 
component of UHECRs; the heavy component, 
assumed to be iron, is identified with the galactic 
component. 

Having identified the extragalactic component 
of the UHECR spectrum, we fit this component 
using the assumption that the protons have lost 
energy between their acceleration and our obser- 
vation, and that the sources are distributed uni- 
formly throughout the universe, with a density 
that is smoothly varying in time. The fits assume 
a simple power law spectrum, with a spec- 

tral index 7 identical for all sources. We use the 
average energy loss model of Berezinsky et aL[S] 



to find the observed spectrum from a source a 
given redshift. 

In addition, we use a Monte Carlo (MC) 
simulationjHI for the discrete energy losses asso- 
ciated with pion production. The MC is only im- 
portant for propagating particles from small red- 
shifts, z < 0.1, where the proton interacts only 
a few times. For propagation from larger red- 
shifts, the distribution of observed energies for a 
given input energy is a gaussian of ever narrowing 
width. 

The MC is run to create an Ein-Eo^t matrix 
with two million entries for propagation from each 
redshift, z, to our observation. The Ei„_-Eout ma- 
trix is calculated every 0.01 step in z for z = 
to 2; = 4. Smaller steps, as small as 0.0001, are 
used below 0.05. Sources are assumed to be uni- 
formly distributed at any given redshift, and with 
a density with which varies as (1 -I- z)™. The 
MC results at individual redshifts are added to- 
gether to form an -Ein-^'out matrix for uniformly 
distributed sources evolving according given value 
of m. These z-integrated matrices were formed 
for m =0 to m =5 in steps of 0.5. 

The fits are performed by varying the nor- 
malization, the spectral slope {7), and the evo- 
lution parameter (to). The goodness-of-fit is 
evaluated using the binned-maximum-likelihood 
method. The comparison spectra for fitting are 
calculated by multiplying a vector containing the 
input spectrum by the Ein-Eout matrix using ma- 
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Figure 1. The composition measurements used in 
the toy model of galactic and extragalactic com- 
ponents in the the UHECR spectrum. HiRes Pro- 
totype/MIA measurements are show as triangles, 
HiRes Stereo measurements as circles. The lines 
above and below the measurements are the results 
of Corsika/QGSJetPUl calculations for protons 
and iron primaries respectively. 



trix multiplication. The input spectrum vectors 
are calculated every 0.05 in the spectral index 7. 
The comparison spectra for intermediate values 
of 7 and m are found by interpolation. 

The result of the fit to the most recent HiRes- 
I[7j and HiRes-II 8 monocular spectra is 7 = 
-2.38 ± 0.04 and m = 2.8 ± 0.3. This fit is 
shown in Figure |21 The quoted uncertainties in- 
clude the effect of the correlation between 7 and 
m. This correlation is shown in Figure |3 Not in- 
cluded in this fit are any systematic uncertainties. 
The largest systematic uncertainty is expected to 
be due to the details of composition used in the 
aperture calculation for the HiRes-II spectrum. 
Changes in the composition tend to move the best 
fit point along the trough of Figure 13 changing 
the m vs. 7 relationship little. 



Figure 2. The HiRes-I and HiRes-II monocular 
spectra, with the result of the best fit spectrum. 
HiRes-I points are shown as red squares, HiRes-II 
points as black circles. The one a upper limit on 
the flux in two more bins for each measurement 
are shown above the highest energy actually ob- 
served. The black fit line is the sum of the galac- 
tic (green) and extragalactic (red) components. 
The parameters of the extragalactic component 
are 7 — —2.38 and m — 2.6. There are 42 data 
points, giving 39 degrees of freedom in the fit. 



2. HiRes Monocular Composition Mea- 
surement 

Because the uncertainty in composition is a 
large part of the uncertainty in the parameters 
in the fit to the spectrum, we were motivated to 
try to improve on the existing composition mea- 
surements, especially below 10^^ eV. To go to low 
energies with HiRes requires using the monocular 
data sets, and we chose the HiRes-II monocular 
data set for this analysis. 

Unfortunately, the limited elevation coverage 
of the HiRes-II aperture biases our Xmax accep- 
tance, with the bias increasing at lower energies. 



3 




2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 

^ -y. Spectral Slope 



Figure 3. A scan of (^21ogL) in m-j space. 
The first four a contours are shown. 



The bias stems from the requirement that we find 
^max within that extent of the shower observed 
in the detector. Events that are closer to the de- 
tector are more hkely to have Xmax above the 
visible range, and thus be cut. Furthermore, be- 
cause lower energies can only be observed close 
to the detector, events at these energies will have 
a larger acceptance bias than those at higher en- 
ergies. This bias precludes performing an elonga- 
tion rate analysis at energies below 10^* cV. 

Instead, we have chosen to fit the X^^ax distri- 
bution in energy bins to a combination of MC 
generated X^ax distributions from proton pri- 
maries and from iron primaries. The Xmax distri- 
bution in each energy bin is stored in a histogram, 
and the data histogram is then compared to a 
mixture of the proton and iron MC histograms 
to find the best fit proportion of each. The fit is 
performed using the binned maximum likelihood 
technique as implemented in the HBOOKjS] rou- 
tine HMCLNL. The uncertainty of the fit in each 
energy bin is taken from the width of — 21ogL as 
fit to a quadratic in the region about the maxi- 



mum. An example of this procedure for the bin 
at 10^^ eV is shown in Figure 0] The result of 
these fits for the data at all energies is shown in 
Figure El 
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Figure 4. The results of a fit to X,„ax at an en- 
ergy of 10^* eV. The top panel shows —2 log L as a 
function of the proton fraction in 1% steps (black 
histogram), and a polynomial fit to the same (red 
line). The bottom panel shows the X^^x distri- 
bution that results from the fit, with the data 
shown as filled circles with error bars, the proton 
distribution in blue, the iron distribution in red 
and the sum in black 



This fit gives only the proportion of proton 
and iron in the final, accepted sample; as such, 
it is not intrinsically interesting, but must be 
corrected to find the proton-iron ratio of UHE- 
CRs entering the atmosphere. However, one can 
test this procedure on an independent MC sam- 
ple where one knows a priori which events are 
protons and which are iron. The result of this 
analysis on the MC sample used to calculate the 
aperture used in the spectrum analysis is shown 
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Figure 5. The raw composition measured in the 
data. 



in Figure The fits agree very weU with what is 
actually in the sample. 

The raw composition measurements are cor- 
rected using the relative acceptances in the pro- 
ton and iron MC samples. These are the same 
MC samples used to provide the comparison 
Xmax distributions used in the fit. The fit proton- 
iron proportion is interpreted as given numbers 
of proton events and iron events. These num- 
bers then are separately corrected for the rela- 
tive acceptance appropriate to each. Finally, the 
proton-iron ratio is recalculated using the cor- 
rected numbers of events of each type. 

For a MC sample, we can again compare the 
corrected proton-iron proportion to the input to 
the MC simulation. This is shown in Figure 
In this case, the input is just the HiRes Proto- 
type/MIA and HiRes Stereo composition mea- 
surements, where the mean X^ax of the data 
is interpreted as representing a proton-iron pro- 
portion given by the distance between the Cor- 
sika/QGSJET proton and iron lines. The input 
is in fact given by linear fits to the mean Xyna.x 



Figure 6. The raw composition measured in the 
MC sample used in the spectrum calculation. The 
black, filled circles are the result of the fits. The 
blue, open circles are the actual fractions using 
the known, a priori information available in the 
MC sample. 



above and below lO^"^ ®^ eV. Again, the corrected 
results of the fits agree very well with the inputs 
to the MC. 

Finally, we show as a very preliminary result, 
the corrected proton-iron ratio for the data. See 
Figure |H1 This result shows a very light compo- 
sition (90% protons) above 10^^'^ eV. Below this 
energy we measure a sharply lower composition, 
but with large uncertainties because of few events 
and a large acceptance correction. The measured 
composition at high energies is in agreement with 
the HiRes Stereo measurement, but closer to be- 
ing all protons. 

The data used in the fits comes from the same 
data set used in the HiRes-H spectrum analysis 
which was fit above. A considerable amount of 
data, about a factor of three more, was collected 
after this data set. This later data was collected 
with lower light thresholds, and as such, will have 
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Figure 7. The corrected composition measured in 
the MC sample used in the spectrum calculation. 
The black, filled circles are the corrected result of 
the fits. The green, open triangles and squares 
are the HiRes Prototype/Ml A and HiRes Stereo 
composition measurements which form the input 
to the MC. 
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Figure 8. The preliminary, corrected composition 
measured in the data. The black, filled circles are 
the corrected result of the fits. For comparison, 
the HiRes Prototype/MIA and HiRes Stereo com- 
position measurements are shown as green, open 
triangles and squares, respectively. 



a lower energy threshold. In addition, the aper- 
ture at each energy will be increased, reducing 
the size of Xmax biases. Performing this analysis 
on this expanded data set should give us much 
tighter constraints on the composition in the mid- 
dle of the 10^^ eV energy decade. 

The Telescope Array (TA), which has already 
received Japanese funding, and the TA Low En- 
ergy extension (TALE) will be able to push 
this analysis even further. TALE will use the 
HiRes phototubes with larger mirrors to push the 
thresholds even lower. The mirrors will be ar- 
ranged in a tower rather than a ring to avoid bi- 
ases in Xmax- Geometrical resolution at the level 
of HiRes Stereo will be provided by coincidences 
with the TA ground array. 
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